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What is claimed is: 



1 . A method for reducing computation time of an analysis of diffraction oi'rncider 
electromagnetic radiation frobi a periodic grating having a direction of periodicity, said 
analysis involving a division if said periodic grating into layers, with an initial layer 
corresponding to a space abovl said periodic grating, a final layer corresponding to a 
substrate below said periodic grating, and said periodic features of said periodic grating 
lying in intermediate layers between said initial layer and said final layer, a cross-section 
of said periodic features being d^scretized into a plurality of stacked rectangular sections, 
within each of said layers a pern^tfmryahd electromagnetic fields being formulated as a 
sum of harmonic componentsialolag saiddirgcjta^h of periodicity, application of 
Maxwell's equations providfag anlintra-layer matrix equation in each of said intermediate 
layers equating a product of a wavl-vector matrix and first harmonic amplitudes of one of 
said electromagnetic fields to a^GOM partial derivative of said first harmonic amplitudes 
of said one of said electromagnetic fields with respect to a direction perpendicular to a 
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plane of said periodic grating, said wave-vector matrix being dependent on intra-layer 
parameters and incident-radiation parameters, a homogeneous solution of said intra-layer 
matrix equation being an expansion of said first harmonic amplitudes of said one of said 
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electromagnetic fields into first exponen|ial functions dependent on eigenvectors and 
eigenvalues of said wave-vector matrix, comprising the steps of: 

\ 

determination of a layer-property parameter region and a layer-property 
parameter-region sampling; \ 
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determination of a maximum harmonic order for said harmonic components 
of said electromagnetic fields; 

calculation of required permitt: vity harmonics for each layer-property value 
in said layer-property parameter region deteifnined by said layer-property parameter- 
region sampling; 

determination of an incident-radiation parameter region and an incident- 
radiation parameter-region sampling; 

calculation of said way^vector biatrix\ased*Qn said required permittivity 

1 

harmonics for said each layer-progferty value 1^ saicjjla^er-propertyjjarameter region 

determined by said layer-propert| parameter-region sampling and for each incident- 

V 

radiation value in said incident-rattiation parameter region determined by said incident- 
radiation parameter-region sampling; 

calculation of eigenvectors and eigenvalues of each of said wave-vector 

matrices for said each layer-property value in saicl layer-property parameter region 

t 

determined by said layer-property parameter-region sampling and for said each incident- 
radiation value in said incident-radiation parameter region determined by said incident- 
radiation parameter-region sampling; 

caching of said eigenvectors and saiid eigenvalues of said each of said 
wave-vector matrices in a memory; and 



use of said eigenvectors and said eigenvalues for said analysis of said 



diffraction of said incident electromagnetic radiation from said periodic grating. 

J 
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2. The method of claim 1 furthoftcomprising the step of caching in, said memory, 
said wave-vector matrices for said eaah layer-property value in said layer-property 
parameter region determined by said layer-property parameter-region sampling and for 
said each incident-radiation value in s lid incident-radiation parameter region determined 
by said incident-radiation parameter-region sampling. 

3. The method of claim 2 further comprising the step of caching in, said memory, 
said required permittivity harmonics for said each layer-property value in said layer- 
property parameter region deterafmed by said\ayer-property parameter-region sampling. 



4. The method of claim/l further cqmprisir 



g the step o£.calculating a product of a 



T 

SS 1 0 square root of each of said Eigenvalues and a corresponding one of said eigenvectors for 

l/p 



said each layer-property valVe in said laW-property parameter region determined by said 



€| layer-property parameter-regioh^sa^fing and for said each incident-radiation value in 
"L said incident-radiation parameter region determined by said incident-radiation parameter- 
ry region sampling. 

15 5. The method of claim 4 further comprising the step of caching in said memory, said 

product of said square root of said each of said eigenvalues and said corresponding one of 
said eigenvectors for said each layer-property value in said layer-property parameter 
region determined by said layer-property p Jameter-region sampling and for said each 
incident-radiation value in said incident-radiation parameter region determined by said 
20 incident-radiation parameter-region sampling! 

6. The method of claim 1 wherein another of said electromagnetic fields is 
expressible as an expansion of second harmoni<^amplitudes into second exponential 
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functions dependent on said eigenvectors and said eigenvalues of said wave-vector 
matrix, application of boundary conditions of said electromagnetic fields at boundaries 
between said layers provides a boundaii-matched system matrix equation, and solution 
of said boundary-matched system matrix equation provides said diffraction of said 
incident electromagnetic radiation from laid periodic grating, and wherein said use of 
said eigenvectors and said eigenvalues f<Sr said analysis of said diffraction of said 
incident electromagnetic radiation from said periodic grating comprises the step of: 

discretization of a crfTss-section of^a ridge of said periodic grating into a 



stacked set of rectangles on said substrate; 



/ 

id .sub: 



■y jo retrieval, from said memory, for each of said rectangles, of said 

Si V / 

O eigenvectors and said eigenvalues\^^dJ^e-vector matrix based on said intra-layer 

?yy j 

parameter values of said each of said rectangles, and based on said incident-radiation 

N I 

^ parameter values of said incident electromagnetic radiation; 

° 

nJ construction of said boundaxt-matched system matrix equation using said 

m 

1 5 eigenvectors and said eigenvalues of said wave-vector matrices retrieved from said 



memory for said each of said rectangles; and] 

solution of said boundary-matched system matrix equation to provide said 
diffraction of said incident electromagnetic radiation from said periodic grating. 

7. The method of claim 1 wherein said intra-layer parameters for one of said layers 

\ 

20 include an index of refraction of a material of said periodic features in said one of said 
layers, an index of refraction of said initial layerAa length of periodicity of said periodic 
features, a width of said periodic features in said oVe of said layers, and an offset distance 
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of said periodic features in said one oi 
include an angle of incidence of said 
electromagnetic radiation. 



said layers, and said incident-radiation parameters 
e&ectromagnetic radiation and a wavelength of said 



8. The method of claim 1 whereh l within said each of said layers, any line directed 
normal to said periodic grating passes through a single material. 



9. The method of claim 1 whereir 
mathematically approximated as semi 

10. The method of claim 1 




said initial layer and said final layer are 
infinite. 

layer-property parameter region and said 



incident-radiation parameter region describe a^yper-rectarigle. 

11. The method of claim 1 whereh coefcclents of said expansion of said harmonic 
amplitudes of said electronkgnetic field into saidexponSTff^fonctions include factors 
which are elements of an eigfeivector matrix obtained from said wave-vector matrix, and 
exponents of said expansion of said harmonic amplitudes of said electromagnetic field 
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include factors which are square roots of eigenvalues of said wave-vector matrix. 

12. The method of claim 1 1 wherein said layer-property parameter-region sampling 
is at a uniform density. 

13. The method of claim 1 1 wherein|said layer-property parameter-region sampling 
is at a non-uniform density. 

14. The method of claim 12 wherein ^aid layer-property parameter-region sampling 
is done on a uniform grid. 

15. The method of claim 12 wherein s^d layer-property parameter-region sampling 
is done on a non-uniform grid. 
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16. The method of claim 1 1 wlkrein at least one dimension of said incident-radiation 
parameter region has a range of a sii igle value. 

17. The method of claim 1 1 whelrein at least one dimension of said layer-property 
parameter region has a range of a single value. 

1 8. A method of determining din^ensions of a physical profile of a 

repeating, regularly-spaced series of sfructures, comprising the steps of: 

illuminating a test area of s^^pre^of said structures with incident 
radiation having a plurality of Wavelengths, s^d incident radiation 
traveling along an optical path; 

measuring radiation diffracted from said test area at said plurality of 
wavelengths to obtain a measured diffraction spectrum; 

determining an intra-l^fe^agAd an intra-layer^ampling of 



intra-layer parameters corresponding to lasers of profile shapes to be included in a first 

I 

sub-library; ii 

determining an incident-radiation range and an incident-radiation 



sampling of incident-radiation parameters corresponding to said optical 
path of said incident radiation; 

generating pre-calculated intra-layer-dependent portions of a diffraction 
calculation for said intra-layer sampling of said intra-layer parameters within said nitra- 



layer range, and said incident-radiation sampling of said incident-radiation parameters 
within said incident-radiation range; 

caching said pre-calculated intra-layer-dependent portions of said 



52 



diffraction calculation for said intra-lapr sampling of said intra-layer parameters within 
said intra-layer range, and said incident-radiation sampling of said incident-radiation 
parameters within said incident-radiatibn range; 

calculating, utilizing said pre-calculated intra-layer-dependent 

5 portions of said diffraction calculation, For each of said profile shapes in said first sub- 
library, a calculated diffraction spectrum according to said diffraction calculation to 
provide a second sub-library of calculated diffracted spectra; 

indexing said first sub-library oO^d^T^file shapes with 
corresponding spectra from said^econd s^ib-Iibraf^ of said calculated diffraction spectre 

1 0 to provide a library of profile-Spectra pairs; anc 



comparing said measpred diffraction speqtrW with said calculated 
diffraction spectra in said libkry of said nppfile-spectra pairs to find a 



best-match calculated spectnim wttehmost closely matches said measured 

\ 

diffraction spectrum, whereby one of said profile shapes corresponding to 



1 5 said best-match calculated spectrum is a best|match for said physical 

profile, and said dimensions of said physical profile are determined to be 

corresponding dimensions of said one of said profile shapes. 

I, 

19. The method of claim 18 wherein said optical path includes a focusing 
mechanism. 

20 20. The method of claim 1 8 wherein said measuring of said diffracted light 

involves a measurement of intensity of said diffracted light. 

1 

21. The method of claim 18 wherein said measuring of said diffracted light 
involves a measurement of intensity of said diff|j*cted light, and a phase 
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difference between polarizations of slid diffracted light. 

22. The method of claim 18 wherein said comparing of said measured 
diffraction spectrum with said calculated diffraction spectra in said 
library to determine a best-match calculated spectrum is performed using a 
least squares computation. 

23. The method of claim 22 wherein said least square computation is 
performed by selecting a number c^jdmpte^avelengths, subtracting said 
measured diffraction spectiWwith said calculated diffraction spectra at 

/ . ...... 

said sample wavelengths tft determine di|jerenc^s at eacn oi saia sample 
0 10 wavelengths, and squaring and summing^^ 

24. An apparatus tor determining dimensions of a physical profile of a 

repeating, regularly-spacecT^rie|^f structures, comprising: 

\ 

an electromagnetic radiation sjource for generating electromagnetic 
radiation having a plurality of wavelengths; 

a focusing means for directing (said electromagnetic radiation along an 
optical path to illuminate a test area of tsaid series of said structures with 
incident radiation; 

a photometer for measuring diffraction of said incident radiation from 
said test area at said plurality of wavelengths to obtain a measured 
20 diffraction spectrum; 

a first memory; 
a look-up table; 

means for generating a first sub-liljrary of profile shapes; 
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means for storing said first s^b-library of said profile shapes in a 
first array of said look-up table; 

means for generating pre-calctalated intra-layer-dependent portions of a 
diffraction calculation for an intra-laMer sampling of 
intra-layer parameters within an intra|ayer range, and an 
incident-radiation sampling of incident-radiation parameters within an 



incident-radiation range, and caching said pre-calculated 



intra-layer-dependent portions of said 



tion calculation in 



said tirst memory; 

means for calculating, iftilizing said prd-calculated intra-layer-dependent portions 

/ \ A 

of said diffraction calculation^ for each of saitf profile shapes in said first 

according to said rigorous 
coupled-wave calculation to pro\^id& I a^f cond sub-libr^y of calculated diffraction 
spectra, 



sub-library, a calculated diffre^tion spectrum 



means for storing said second sub '^library of said calculated diffraction spectra in 
a second array of said look-up table at entries corresponding to corresponding ones of 



said profile shapes in said first sub-library;^ and 

means for comparing said measured diffraction spectrum with said 



calculated diffraction spectra in said second array of said look-up table to 

find a best-match calculated spectrum which most closely matches said 

v 

measured diffraction spectrum, whereby one of said profile shapes 
corresponding to said best-match calculated spectrum is a best match for 
said physical profile, and said dimensions of said physical profile are 
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determined to be corresponding dimensions of said one of said profile 
shapes. 

25. A method of determining Simensions of a physical profile of a 



repeating, regularly-spaced series of s= 
illuminating a test area of said 
radiation having a plurality of angles a 
traveling along an optical path; 



ctures, comprising the steps of: 
series of said structures with incident 
a single wavelength, said incident radiation 




measuring radiation diffrac^. 
angle to obtain a measured diffraction sSpectrum. 
determining an intra-layer range 



test area at said plurality of 



intra-layer parameters corresp 
sub-library; 




and ^ inif a-layer sampling of 

Tofile shapes to be included in a first 



determining an incident-radiationjrange and an incident-radiation 
sampling of incident-radiation parameters; corresponding to said optical 



t 



path of said incident radiation; 

generating pre-calculated intra-layer-dependent portions of a rigorous 
coupled wave calculation for said intra-layer sampling of said intra-layer 
parameters within said intra-layer range, and said incident-radiation 
sampling of said incident-radiation parameiers within said 
incident-radiation range; 

■I 

caching said pre-calculated intra-layer-dependent portions of said 
diffraction calculation for said intra-layer sampling of said 
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intra-layer parameters within said intA-layer range, and said incident-radiation sampling 
of said incident-radiation parameters vdthin said incident-radiation range; 

calculating, utilizing said pre-cdculated intra-layer-dependent 
portions of said diffraction calculation for each of said profile shapes in said first sub- 
library, a calculated diffraction spectrun according to said diffraction calculation to 
provide a second sub-library of calcula ed diffracted spectra; 

indexing said first sub-library of said profile shapes with 
corresponding spectra from said second sub-library of said calculated diffraction spectra 
to provide a library of profile-spectra p£ irs; 

comparing said measured diffraction sM^trum with said calculated 
diffraction spectra in said liprary of said|pro|fli-sj>ectra pairs to find a 
best-match calculated spectkm which most closehwngtches said measured 
diffraction spectrum, wherebybBS-tff said profile shapes corresponding to 
said best-match calculated spectrum is a b'est match for said physical 



profile, and said dimensions of said physical profile are determined to be 



0^1 



corresponding dimensions of said one of said profile shapes. 



I 



26. The method of claim 25 wherein said optical path includes a focusing 



mechanism. 1 

27. The method of claim 25 whereinlsaid measuring of said diffracted light 
20 involves a measurement of intensity of said diffracted light. 

28. The method of claim 25 wherein said measuring of said diffracted light 
involves a measurement of intensity of said diffracted light, and a phase 
difference between polarizations of said diffrabted light. 
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29. The method of claim 28\wherein said comparing of said measured 
diffraction spectrum with said calculated diffraction spectra in said 
library to determine a best-match calculated spectrum is performed using a 
least squares computation. I 

30. A method for analyzing tHe profile of a workpiece comprising: 
generating intra-layer-dependent portions of diffraction calculations for a plurality 

of profile shapes; 

caching said intra-layer-depenAent portions of said diffraction calculations in a 
memory; 

10 generating calculated-dif|factio| spectr^for said profile shapes based on said 

intra-layer-dependent portions j£>f said d^ffractioi^calculations retrieved from said 
memory; 

generating a library of ssftd profile! shapes anchsqid calculated-diffraction spectra; 

illuminating a portion of the workpiece with incident radiation; 
1 5 measuring radiation diffracted froin said portion of the workpiece; 

obtaining a measured-diffraction spectrum based on the measurement of the 
diffraction radiation; 

matching said measured-diffraction spectrum with one of said plurality of 
calculated-diffraction spectra contained in said library; and 
20 selecting the profile shape associated with the calculated-diffraction spectrum that 

matches said measured-diffraction spectrum. 

3 1 . The method of claim 30 further comprising intra-layer parameters associated 
with said profile shapes and incident-radiation parameters associated with said incident 
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radiation, and wherein said step of generating intra-layer-dependent portions further 



compnses: 

determining an intra-layer range and an intra-layer increment of said intra-layer 
parameters; 

determining an incident-radiation fange and an incident-radiation increment of 
said incident-radiation parameters; and 

generating said intra-layer-dependept portions of said diffraction calculations for 
said intra-layer parameters between said in xa-layer range at said intra-layer increment, 
and for said incident-radiation para#€ters qetwedq said incident-radiation range at said 
1 0 incident-radiation increment. 

32 The method of claim 30 wherehi said|i^asuring step further comprises 

| 1 . 

measuring the intensity of said 'diffracted radiafroi 



15 




33. The method of claim\2 furth^omprism^ the step of measuring a phase 
difference between said incident radiation and said diffracted radiation. 



34. The method of claim 30 wherein said matching step is performed utilizing a 



least square computation. 



35. The method of claim 34 further comprising: 
selecting a plurality of sample wavelengths; 

obtaining a difference at each sample wavelength by subtracting said measured- 
20 diffraction spectrum with said calculated-diffraction spectrum at said sample wavelength; 
and 

squaring said difference; and 
summing said difference. 
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36. The method of claim 30 v herein said illuminating step further comprises 
focusing said incident radiation. 

37. The method of claim 36 wnerein said incident radiation includes light having 
a plurality of wavelengths, and whereil the portion of the workpiece is illuminated from 
a single angle. | 

38 The method of claim 3 0 wherein said incident radiation includes light having 

t 

a single wavelength, and wherein the portion of the workpiece is illuminated from a 

plurality of angles. ; 

39. A method of analyzing t^^r^oiog^f a semi-conductor wafer utilizing an 
incident light having a plurality of^wavelengths, said method comprising: 

providing a library of a plurality of profile dmpes and a plurality of calculated- 

/ \ n \ 

diffraction spectra, wherein each profile shape is ^fjso^iated with each calculated- 
diffraction spectrum; 

focusing the incident light; 



illuminating a portion of toe^wafer from a single incident angle; 

\ 

measuring light diffracted from said portion of the wafer; 



obtaining a measured-diffraction spectrum based on the measurement of the 
diffracted light; 

comparing said measured-diffraction spectrum with said calculated-diffraction 

spectrum; and | f 

obtaining a profile for the wafer based onWd comparison. 
40. The method of claim 39 wherein said generating step further comprises: 
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determining an intra-layer ran^ ; and an intra-layer increment of intra-layer 



parameters associated with said profile 



shapes; 



determining an incident-light range and an incident-light increment of incident- 
light parameters associated with said ir.cident light; 

generating intra-layer-depender.t portions of diffraction calculations for said intra- 
layer parameters between said intra-layer range at said intra-layer increment, and for said 
incident-light parameters between said incident-light range at said incident-light 
increment; and 

generating said calculated-diffraction specixa based on said intra-layer dependent 




comparing step further comprises: 



portions. 

41. The method of els 
selecting a plurality of \ 

obtaining a difference at ea<^isajjpfffe wave^ngth by subtracting said measured- 
diffraction spectrum with said calculated-jdiffraction spectrum at said sample wavelength; 
and 

squaring said difference; and 
summing said difference. 

42. The method of claim 39 wherein said measuring step further comprises: 
measuring the intensity of the diffracted light; and 

measuring a phase difference between said incident light and said diffracted light. 

43. A method of analyzing the prof^e of a workpiece utilizing an incident light 
having a single wavelength, said method comprising: 

generating a library of profile shapes;\ 
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determining an intra-layer range and an intra-layer increment of intra-layer 
parameters associated with said profile siapes; 

determining an incident-light range and an incident-light increment of incident- 
light parameters associated with said incident light; 
5 generating intra-layer-dependent portions of diffraction calculations for said intra- 

layer parameters between said intra-layer range at said intra-layer increment, and for said 
incident-light parameters between said incident-light range at said incident-light 
increment; 

generating a plurality of calcui^fed-pffractior^ spectrum based on said intra-layer 
O 10 dependent portions; 

M storing said plurality of calcplated-diffractioy plectrum in said library of profile 

shapes, wherein each calculated-dif^action spectpdm \s]associated with each profile 
shape; 

illuminating a portion of the wafer from a plurality of incident angles; 

1 

measuring light diffracted from said portion of the wafer; 
obtaining a measured-diffraction spectrum based on the measurement of the 
diffracted light; 

comparing said measured-diffraction spectrum with said calculated-diffraction 
spectrum; and 

20 obtaining a profile for the wafer based on said comparison. 

44. The method of claim 43 further com|rising: 

caching said pre-calculated intra-layer-de|endent portions of said diffraction 
calculations. 
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45. A method of analyzinathe profile of a workpiece utilizing an incident light, 
wherein the workpiece has a plurality of layers, said method comprising: 
generating a library of profile shapes; 

determining an intra-layer rarlge and an intra-layer increment of intra-layer 
5 parameters associated with said profile shapes; 

determining an incident-light range and an incident-light increment of incident- 
light parameters associated with said incident light; 

generating intra-layer-dependenlportions of diffraction calculations for said intra- 
layer parameters between said intra-lay^crange afraid intra r layer increment, and for said 
1 0 incident-light parameters between s&d inlident-light Wge at said incident-light 
increment; 

generating a plurality of ealculated-diffractioh Spectrum based on said intra-layer 
dependent portions; 

storing said plurality of calculale^^ffraction spectrum in said library of profile 
15 shapes, wherein each calculated-diffraction spectrum is associated with each profile 

shape; \ 

illuminating a portion of the workpiece v^th said incident light; 
measuring light diffracted from said portio^L of the workpiece; 
obtaining a measured-diffraction spectrum based on the measurement of the 

20 diffracted light; ^ ^ 

comparing said measured-diffraction spectrum with said calculated-diffraction 

spectrum; and \ 

obtaining a profile for the layers of the workpiece based on said comparison. 
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46. The method of claim 45 wherein said comparing step further comprises: 
selecting a plurality of sample wa\felengths; 

obtaining a difference at eadh sample wavelength by subtracting said measured- 
diffraction spectrum with said calculated-diffraction spectrum at said sample wavelength; 
and 

squaring said difference; and 
summing said difference. 

47. The method of claim 45 whepettTSaid measuring step further comprises: 
measuring the intensity of Aje diffracted light; and . ( 

measuring a phase difference between said inbident light and said diffracted light. 

r orkpiece comprising: 



;d to illuminate a portion of the 



48. A system for analysing the profile of a 
an electromagnetic radiation source confi^ur 
workpiece with an incident radikion; 

a photometer configured to measure radiation diffracted from the workpiece; 

\ 

a library having a plurality of profile shapes and a plurality of calculated- 

\ 

diffraction spectra, wherein said calculated-diffraction spectra were generated based on 

\ 

intra-layer dependent portions of diffraction calculations of said profile shapes; and 



a processor configured to compare said measured-diffraction spectrum with said 

\ 

calculated-diffraction spectrum to obtain a profile qf the workpiece. 



49. The system of claim 48 further comprising a cache of said intra-layer 

\ 

dependent portions of diffraction calculations, wherein said intra-layer dependent 
portions were generated for intra-layer parameters associated with said profile shapes, 
between an intra-layer range at an intra-layer increment* and for incident-light parameters 
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associated with said incifleVt light, between an incident-light range at an incident-light 
increment. 

50. The system of cla\m 48 further comprising a focus lens configured to focus 
said incident radiation. 

5 51. The system of claim^ wherein said incident radiation includes light having 

a plurality of wavelengths. 

52. The system of claim 48 wherein said incident radiation includes light having 
a single wavelength. 

53. The system of claim 48 whlbdn saichjAotometei; is configured to measure the 
10 intensity of said diffracted radiation/ 

S / 

J| 54. The system of claim 58 wherein said photometer is configured to measure a 

2 phase difference between said incident radiation adjd said diffracted radiation. 

%\ 55. The system of claim 48 wherein said processoTisfifftKer configured to 

□ perform a least square computation. 

U 15 56 A system for analyzing the profile of a workpiece utilizing an incident light, 

5 \ 

i wherein the workpiece has a plurality of layers, said system comprising: 
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means for generating a library of profile shapes; 

means for determining an intra-layer range and an intra-layer increment of intra- 
layer parameters associated with said profile shapes; 
20 means for determining an incident-light range and an incident-light increment of 

incident-light parameters associated with said incident light; 



means for generating intra-layer-dependent portions of diffraction calculations for 
said intra-layer parameters between said intra-layer range arsaid intra-layer increment, 
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and for said incident-light parameters between said f cident-hght range at said incident- 
light increment; 

means for generating a plurality of calculated-diffraction spectrum based on said 

intra-layer dependent portions; 

means for storing said plurality of calculated-diffraction spectrum in said library 
of profile shapes, wherein each calculated-djffr^n spectrum is associated with each 
profile shape; 

means for illuminating a portion of m workpiec^ with said incident light; 
means for measuring light aiffracled frpm said portion of the worlcpiece; 
means for obtaining a measure£atfpaction spectrum based on the measurement 

of the diffracted light; 

means for comparing said measured-diffraction spectrum with said calculated- 
diffraction spectrum; and 

means for obtaining a profile for /he layers of the workpiece based on said 



comparison. 
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